On the Possibility of Measuring the Abraham Force using Whispering Gallery Modes 
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Critical experimental tests of the time-dependent Abraham force in phenomenological electrody- 
namics are scarce. In this paper we analyze the possibility of making use of intensity-modulated 
whispering gallery modes in a microresonator for this purpose. Systems of this kind appear at- 
tractive, as the strong concentration of electromagnetic fields near the rim of the resonator serves 
to enhance the Abraham torque exerted by the field. We analyze mainly spherical resonators, al- 
though as an introductory step we consider also the cylinder geometry. The order of magnitude 
of the Abraham torques are estimated by inserting reasonable and common values for the various 
input parameters. As expected, the predicted torques turn out to be very small, although probably 
not beyond any reach experimentally. Our main idea is essentially a generalization of the method 
used by G. B. Walker et al. [Can. J. Phys. 53, 2577 (1975)] for low-frequency fields, to the optical 
case. 

PACS numbers: 42.50.Wk, 42.50.Tx, 03.50.-z 
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I. INTRODUCTION 

The one-hundred years old Abraham-Minkowski 
energy-momentum problem in phenomenological electro- 
dynamics P, [1] has recently attracted considerable inter- 
est. Assume henceforth for simplicity that the medium is 
nonmagnetic and nondispersivc, with refractive index n. 
In our opinion - as expressed in the review article some 
years ago by one of the present authors f3] - the most 
physical expression for the electromagnetic force density 
is the Abraham expression (SI units assumed) 



pAM 



1 d 



dt 



(E X H) 



(1.1) 



Here the first term f^'^ = -{eo/2)E'^Vn'^ is different 
from zero in regions where n varies with position, espe- 
cially in the surface regions of dielectrics. This term is 
common for the Abraham and Minkowski tensors, and 
may appropriately be called the Abraham-Minkowski 
term. The second, time-dependent term in Eq. (jl.ll) . is 
the Abraham term. It may be noted that the expression 
(II. ip is in agreement with Ginzburg as well as with 
Landau and Lifshitz [5]. 

One may ask: is it possible to detect the Abraham 
term in experiment? The answer is yes, but the task has 
proven to be surprisingly difficult. The magnitude of the 
electromagnetic frequency is a significant factor in this 
context. Let us give a brief account of three important 
experimental cases: 

1) The first case is the quasi-stationary torque exper- 
iment of Walker et al. [1, Q- Strong, time- varying, or- 
thogonal electric and magnetic fields were applied across 
a dielectric shell of high permittivity, making it possi- 
ble to detect the oscillations themselves. In this way the 
Abraham term was measured quantitatively. 
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2) When considering instead high-frequency fields such 
as in optics, the Abraham term fluctuates out when av- 
eraged over a period. One can thus no longer detect this 
force directly. The physical effect of this force is how- 
ever to produce an accompanying mechanical momen- 
tum propagating together with the Abraham momentum. 
The resulting total momentum is the Minkowski mo- 
mentum, corresponding to the divergence- free Minkowski 
energy-momentum tensor. This tensor has the partic- 
ular property of being space-like, corresponding to the 
possibility of getting negative field energy in certain in- 
ertial frames. An authoritative experiment measuring 
the Minkowski momentum is that of Jones et al. 8, 9], 
measuring the radiation pressure on a mirror immersed 
in a dielectric liquid. Both cases 1) and 2) are discussed 
in some detail in Ref. 3;]. 

3) The third example to be mentioned is the photon 
recoil experiment of Campbell et al. jlQ] , where the pho- 
ton momentum in a medium (in this case a Bose-Einstein 
condensate) was found to be equal to the Minkowski 
value ?ik. 

Most other experiments are measuring not the Abra- 
ham term but rather the surface force f^'^, although 
claims are sometimes made to the contrary. In our opin- 
ion this is the case also for the interesting new fiber op- 
tical experiment of She et al. [ll|; cf. the remarks in 
Refs. [l2,[il- 

Our main purpose in the present paper is however not 
to interpret already existing experiments, but instead to 
propose the idea of using whispering gallery modes as a 
convenient experimental tool to detect the Abraham term 
in optics. To our knowledge this idea has not been con- 
sidered before. Whispering gallery modes are commonly 
produced in microspheres; they have a large circulating 
power, about 100 W typically, and the field energy is 
concentrated along the rim of the sphere. That means, if 
such a sphere is suspended in the gravitational field and 
fed with an appropriate intensity modulated field, the 
sphere becomes exposed to a vertical torque according 
to Eq. (|1.1[) . With the field energy essentially concen- 
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trated along the rim, the arm in the torque calculation 
is essentially the same as the radius, thus maximizing 
the torque. In effect, this is the idea of the experiment 
of Walker et al. [1, @|) generalized to optical frequencies. 
We have actually suggested this idea qualitatively before, 
in Refs. [Ulil. 

The next two sections give quantitative estimates for 
performing such an experiment. The torque turns out 
to be small, as expected, but not beyond any possibility 
for experimental detection. Spherical geometry, as men- 
tioned, is most typical for the whispering gallery setup. 
In the next section we consider however as an introduc- 
tory step the somewhat more simple geometry of a cylin- 
drical shell. 

Before closing this section, let us give a few more ref- 
erences to the Abraham-Minkowski problem, in addition 
to the references given above. A nice introduction can 
be found in M0ller's book (TJ]. A review, up to 2007, is 
given by Pfeifer et al. p^ . Some more recent papers are 
Refs. [iB-Ei- 



II. CYLINDRICAL GEOMETRY 

Consider first as the simplest case a compact cylinder 
of length L and radius a. On the inside, r < a, the 
permittivity is e and the permeability /i. On the outside, 
r > a, a vacuum is assumed. The dispersion relation for 
stationary modes is known to be [l9| 



region 11911: 



Jrn (u) 



(2.1) 



We are working with SI units and let e and /i be dimen- 
sional, so that D = eE, B — jiH. The transverse wave 
vectors on the inside and the outside are 



Al — nuj/c, A2 — ll>/c, 



(2.2) 



respectively, while their nondimensional counterparts are 



u = Xia, V = A2a. 



(2.3) 



An important property of this equation is that when 
the axial wave vector fc = - as is of interest here as we 
we consider azimuthal modes only - the right-hand side 
vanishes and the problem becomes separable into TE and 
TM modes. 

We write the mode expansions for the fields in the inner 



m— — 00 



00 

^ E 

m— — 00 



'^m(Air) Om-^m; 



and 



He =— ^ Jm(Air)a,„F„ 

m— — 00 
00 



m— — 00 



where 



F„ 



(2.4a) 
(2.4b) 
(2.4c) 

(2.5a) 
(2.5b) 
(2.5c) 

(2.6) 



The coefficients am and bm, corresponding to the TM 
and TE modes, give the weight of each mode. 

In our considerations below we will for simplicity ex- 
tract one single TE mode of high order m, such that there 
is an azimuthally moving momentum concentrated in the 
vicinity of the boundary r — a. (In reality, the incident 
power may be distributed over a band of neighbouring 
m modes, but this does not infiuence the essence of our 
argument.) We first need to determine the magnitude of 
the radial argument Air w u. Let us take 



m = 100, n — 1.5, a — 100 fim. 



(2.7) 



It is known that for a large value of the order m the first 
maximum of the function Jm{x) occurs when x is very 
close to m. This maximum is the one of interest here. 
Thus the lowest resonance frequency uj is determined by 
the equation 



nauj / c — m. 
With the numbers given above, 

w = 2 X 10^^ s~^ 



(2.8) 



(2.9) 



In this manner we manage to make the beam strongly 
concentrated near the rim, as desired. One has in this 
case i?2 = 0, Hr = 0, while the nonvanishing field com- 
ponents of interest are 



Hz =Jm{Xir)bmFm- 



mJm{Xir) braF„ 



(2.10) 
(2.11) 
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The azimuthal component of the Poynting vector S(r) in 
the interior is 

So{r) = -^mrK] - f^J,l(Air)|6,„p, (2.12) 
corresponding to the azimuthal power 

P = L 



S,dr^^^\b^? r^Jlix). (2.13) 



2\i 



In om' case the factor \/x can be extracted outside the 
integral, so that 



2A?zi ' 



dxJl,{x). 



(2.14) 



Assume now that the beam is intensity modulated with 
a frequency ujq (wq low compared with optical frequen- 
cies), 

P Pq cosoiQi, So — Sq cos ujat. (2-15) 
Then the azimuthal Abraham force density is 



J ch — 



1 



dt 



wqS'o sincjo^, (2-16) 



giving rise to the following Abraham torque around 
the vertical symmetry axis: 



Nt 



2ttL / r^f^dr « 2nLa^ / f^dr. (2.17) 
Jo Jo 



Defining the quantity K as 



K = 



-27ra^Po, 



(2.18) 



we thus see that the torque can be written as 

= KujosinuJot. (2.19) 

As expected, the torque becomes very small. As order of 
magnitude we get 



K ^ ^Po ^ (0.7 X 10-2V) • Po 



and the Abraham torque is estimated as 
- (0.7 X lO^^V) • wqPo. 



(2.20) 



(2.21) 



Insert first the very low value of wq ~ 1 s ^, and take 
Po - 100 W. We get - 0.7 x lO'^^ N m, which 
is much less than the value 10~^^ N m obtained in the 
classic Beth experiment [20|, for example, in which the 
angular momentum of light was measured. It is however 
possible to improve the situation by exploiting the fact 
that the build-up and ringdown times for this kind of 
resonators are known to be very small, in the order of 
tens to hundreds of Ns (see discussion below). It is thus 
realistic to insert a much higher value for loq. Inserting 



tentatively ujq = 1000 s'^ we get ~ 0.7 x 10"^^ N m, 
which is perhaps not so unrealistic after all. 

It is physically instructive to look at the system in 
another way, by considering the angular deflection of 
the cylinder instead of the magnitude of the torque. Let 
the cylinder be hanging vertically in the gravitational 
field, suspended by a thin wire of known torsion constant 
K. Denoting the eigenfrequency of the cylinder in the 
absence of any torque by fi, and denoting the damping 
coefficient by 7, we have as equation of motion 



2 K 
- + fl (p = —loq sinujQt. 



(2.22) 



Here / — ^Ma^ is the moment of inertia about the z axis, 
M — pa?L being the cylinder mass with p the material 
density. In our notation, n = IVl^ . With a — 100 |^m as 
above we obtain, when choosing L = 1 mm and assuming 
p ^ 10'^ kg/m-^. 



(2.23) 



For the magnitude of k we may choose a typical value 
characteristic for torsion experiments testing the equiva- 
lence principle, k ~ 10~^ N m/rad [2l|, [23. Then, 



Q. - 10-^ rad s" 



(2.24) 



The magnitude of f2 is large because a is assumed small. 

The largest oscillations occur at resonance, when ljq is 
chosen equal to VI. Then, 



K 
'1^ 



■ cosQt. 



The maximum value, when Pq ^ 100 W, is 



71^-1 4 Pn 



^max — n 



M 7 



10"^ rad 

7 s 



(2.25) 



(2.26) 



It would be of interest to make an estimate of the damp- 
ing constant 7 here, but we postpone that until the next 
section. 

Notice that the very existence of a oscillatory move- 
ment would be enough to make the experiment critical 
with respect to the Abraham force. The Minkowski ten- 
sor does not predict there to be an azimuthal movement 
at all. 



III. SPHERICAL GEOMETRY 

As mentioned above, whispering gallery modes are usu- 
ally associated with microspheres. Let the radius of the 
sphere be denoted by a. As above, we look for the eigen- 
modes, and we will for simplicity focus on the TE modes 
only. (The meaning of the symbol TE is here that the 
electric field is transverse to the radius vector r.) We 
introduce quantities a and r defined by 



Loa/c, 



r/a. 



(3.1) 
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Thus a is the magnitude of the nondimensional wave 
vector in the exterior region (vacuum), whereas f = 1 at 
the boundary. Making use of the Riccati-Bessel function 



= xji{x), 



(3.2) 



the basic TE modes in the interior can conveniently be 
written as 

Er =0, (3.3a) 
Eg = z -T—;^ — Tpi[nar) Fm, (3.3b) 



Erf, 



nar sm 
Ai^ dP[^ {cos 

nar dO 



ipi (nar) F„ 



(3.3c) 



and 
He 

Ha, 



1(1 + 1) Ai lpm(^pgg)^^(^^-)j,^^ (3 4^) 
1 Ahn ldP™(cos( 



za;/i r a 



dd 



m Aim 1 



■ ij'i{naf)Fm, (3.4b) 



P/" (cos 9) ^j'i{naf)Fm, (3.4c) 



ti;/.tsmf r a 

where Aim are constants, and 

Fm = e™-^-*"*. (3.5) 

The mode expansions above essentially follow Stratton 

The components of Poynting's vector are, when aver- 
aged over an optical period. 



Sr ^-^[EeH^ ~ E,pHg], 
Se ^2^[E,i,H*], 
S,p^~^^[EeH:]. 



(3.6a) 
(3.6b) 
(3.6c) 



Assume that the sphere is fed by an incident flux from 
the outside such that only the component 5*0 of S in the 
interior is different from zero. With an intensity modu- 
lated energy flux such as above, S'^ = Sq coswoi, we thus 
get for the azimuthally directed Abraham force density 
in the interior 



-1 



-ujqSq sin Wot. 



(3.7) 



From the above expressions, 

m /(/ + 1) \Aim\'[Pr? ,2 



V'i- (3.8) 



The Abraham torque, directed along the z axis, then 
becomes 



Jirx f^),dV = J rf^smOdV, (3.9) 



where the integration is over the sphere, with dV = 
sin edrd0d(l). Making use of Eqs. ((XT)) and ((X^ we 
obtain 



— 1 Tnna^ l{l + 1) 



where Ki and Ku are the integrals 



Aim\ ujqKiKu sin Wot, 



(3.10) 



Ki = ipf{naf)df 
Jo 

= [V'f - 'ipi-i{na)'ipi+i{na)] , (3.11a) 



Ku = / [P;™ (cos 6*)] 2 sin 6*^6' 
Jo 

2 (/ + m)! 



2/ -1-1 {I -my: 



(3.11b) 



We want to relate this to the total power P flowing in 
the azimuthal direction in the sphere. We calculate P by 
integrating 5*0 over the area of a semicircle with radius 
a, 



P 



where 



d9 / rdrSijj 

^0 

ma + 2 



2{naY ojfjL 



AimVKuiKiY cos unt, (3.12) 



Km = [ ^^f{naf) 
Jo ^ 



Kw = 



[P™(cos^ 
sin 9 



-d9. 



(3.13a) 
(3.13b) 



As before, it is assumed that the supplied power is inten- 
sity modulated, P — Pq coswo^- 

The two last integrals can be processed further, at least 
approximatively. First, we can rewrite Km as 



1 /■"" dx , 
^iii = -nna / — J^{x), 
^ Jo 



(3.14) 



with V = 1 + 1/2. For actual physical values, na ^ 1. We 
can thus replace the upper limit with infinity, and make 
use of formula 6.574.2 in Ref. [23| to get 



nna 



III 



2(2/ + 1) 



(3.15) 



Finally, the integral Kiy is simply (cf. formula 8.14.14 
in Ref. Q) 



IV 



(/ + m)! 
m[l — m)\ 



(3.16) 



We are now able to relate the torque to the power 
P. The result becomes quite simple: 

A n^ — 1 Ama?ujr, 

N^^ 5 ^Posinwoi 

na 

X [i/jfina) - ^i-i{na)iJi+i{na)]. (3.17) 



5 



The radius of the sphere is seen to appear in the prefactor 
a^, as well as in the nondimensional parameter a = uja/c. 
The parameter I occurs only as an order parameter in the 
function V'z- We see that the torque is proportional to m. 
This is as we would expect, as the whispering gallery 
modes are associated with m — I, i.e. the maximum 
value of m. It should correspond to a maximum angular 
momentum and accordingly a maximum torque. 

To proceed quantitatively, the value of a has to be 
determined. For the TE modes it is determined by the 
dispersion relation [l^ 



^^. ipi{na) 



(3.18) 



where £!l^\x) — xh^i^\x) is another member of the 
Riccati-Bessel functions. The equation (|3.18p is complex 
and does not in general have real solutions, but approx- 
imate solutions with only a small imaginary inequality 
are found close io a~ I for ^ ^ 1. 

As at the end of the previous section, we focus now 
attention on the magnitude of the angular deflection 
as this is most likely the quantity of main experimental 
interest. Without changing the notation we write the 
Abraham torque in the form ~ Kojo sin ojot as before, 
where now 



K 



1 Ama? 



(? na 
X [^l'f{na) - il;i^i{na)'iJji+i{na)]Po 



(3.19) 



The equation of motion for takes the same form (j2.22|) 
as before, where now the moment of inertia is 



2 2 r 

5 15 



(3.20) 



M being the mass of the sphere. For definiteness let us 
take a = 100 /im. Then, with p ^ 10"^ kg/m^ we get 
M « 4 /ig and so, with k ^ 10^^ N m/rad as before. 



n - 10** Vk-- 10^ rad s-^ 



(3.21) 



need to know the viscous torque on a sphere executing 
rotary oscillations about its symmetry axis. The solu- 
tion of this problem is shown in Ref. psj . An impor- 
tant parameter in this context is the penetration depth 
S — ■y/2i'/r2, where v is the kinematic viscosity of the sur- 
rounding medium. For air, u = 1.5 x 10~^ m^/s. Thus 
with n ^ 10^ rad s^^ we get S ^ 170 |j.m, which is of 
the same order as a. Strictly speaking we should there- 
fore have to use the complete expression for the viscous 
torque, which is somewhat complicated. For our order- 
of-magnitude considerations it is however sufficient to use 
the simple expression 



(A^.)v 



(3.23) 



(corresponding mathematically to the a/5 1 limit), 
where 77 = 1.8 x lO'^Pa s is the dynamic viscosity for 
air. Identifying (A'^z)viscous with 77 in accordance with 
Eq. (|2.22p . we get for the damping coefficient 



7 = ^ 30 s-\ 



(3.24) 



and the expression p.22p for the maximum deflection can 
finally be written as 



m n — 1 
Zimarja 

X [ipfina) - ^pl-l{na)'^pl+l{na)]Po■ 



(3.25) 



As expected, the deflection is very small. Whereas nu- 
merical evaluation of the ipi functions in general is called 
for, as mentioned, we may note that in cases where I <C 
na the approximation ipi{na) « sin(nQ; — ^) is useful. 
One can moreover obtain a simple estimate of the mag- 
nitude in the cylinder case by inserting 7 from Eq. p.24|) 
into Eq. (I2.26p . whereby one finds 0max ~ 10~* rad. 
Careful adjustments of input parameters are obviously 
needed if the effect is to be verified experimentally. 



IV. ON THE MAGNITUDE OF TORQUES IN 
EXISTING EXPERIMENTS 



With these numerical choices the value of becomes of 
the same order as in the cylinder case. The magnitude 
'?^'max of the maximum deflection at resonance = is 
now 



^ 10m -n? ~1 

Mna 



X -'-2 



■4'i^i{na)ipi+i{na)] — . 

7 



(3.22) 



As we have assumed / ^ 1 and na ^ 1 but otherwise left 
the ratio of these quantities unspecified, the ipi functions 
ought to be calculated numerically. 

Let us finally make an estimate of the magnitude of 
the damping coefficient 7, assuming for definiteness that 
the damping is due to the viscosity of air only. We then 



We close this investigation by making some estimates 
of radiation torques on spheres, as well as on ring res- 
onators (a closely related geometry), for already exist- 
ing experiments. As first example we take the setup re- 
ported in Ref. [1^ , where an infrared laser of wavelength 
A — 1500 nm was used. Two different sphere radii were 
investigated, a = 40 |^m and a = 70 yuon, corresponding 
to values oi a ^ I — m equal to 162 and 283, respectively. 
Although the feeding laser had a power in the order of 
tens of microwatts to milliwatts, the extremely high Q 
factor of the silica sphere meant the buildup of circu- 
lating modes in the sphere grew enormous. Circulating 
powers in excess of lOOW are routinely reported in such 
systems (e.g. [l^l) (although this quantity was not ex- 
plicitly given in the reference (26{). The refractive index 
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of materials used for ultra- high- Q spherical resonators, 
such as fused silica [H, [2^ and quartz [2^, are about 
n = 1.5. With these values as input for Pq we obtain the 
torques [N^ ~ sin uJot] 

f 4 X lO^^'^Nms • Wo r f 40^nl ,,,, 

No~<, in-23AT : for a = <^ . (4.1) 

[ 1 X 10 "^-^Nnis ■ luq [ 70iJ,m ^ ^ 

Note in general that for a sphere, oc a according to 
Eq. p.l7p . whereas 0max oc according to Eq. (I3.25P 
when the viscous dampin g is accounted for. 

The geometry of Ref. [27], which reports circulating 
powers in excess of 100 W, employs the toroidal ring res- 
onator. This geometry has the benefit of having smaller 
mass and therefore smaller moment of inertia than a 
sphere of the same radius, allowing for larger angular 
deflections. For a thin ring the moment of inertia is 

/toroid ~ 27rpAa^ (4.2) 

where A is the area of cross-section. The torque on such 
a toroid would be roughly similar to that on a sphere, 
so it is reasonable to assume the angular deflection to 
be larger and scale as a~^. This could allow larger radii 



which could be beneficial for detection. Cf. also the 
review article [soj . 

We wish finally to re-emphasize the possibility of us- 
ing quite high frequencies loq in order to produce measur- 
able values for the Abraham torque. We assumed above 
the strong field inside the microcavity to react instanta- 
neously to the sinusoidal variations of the input signal, 
an approximation which is good provided the build-up 
and ringdown time (t) of the resonator is small com- 
pared to 2tt/u!o. For the 45|xm radius toroidal resonator 
in Ref. [Blj , for example, a ringdown time of about 43 ns 
was measured. For cavities of even higher Q-factor, ring- 
down times are somewhat longer, yet this implies that we 
may choose tuning frequencies ujq as high as 10^ without 
invalidating the theory. Due to the proportionality of the 
torque with ujq , going close to the megahertz regime could 
increase the torque to perhaps 10~^^ Nm for a sphere 
with radius of some tens of microns. 
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